Spontaneous spiking in an autaptic Hodgkin-Huxley set up 
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The effect of intrinsic channel noise is investigated for the dynamic response of a neuronal cell 
with a delayed feedback loop. The loop is based on the so-called autapse phenomenon in which 
dendrites establish not only connections to neighboring cells but as well to its own axon. The 
biophysical modeling is achieved in terms of a stochastic Hodgkin-Huxley model containing such a 
built in delayed feedback. The fluctuations stem from intrinsic channel noise, being caused by the 
stochastic nature of the gating dynamics of ion channels. The influence of the delayed stimulus is 
systematically analyzed with respect to the coupling parameter and the delay time in terms of the 
interspike interval histograms and the average interspike interval. The delayed feedback manifests 
itself in the occurrence of bursting and a rich multimodal interspike interval distribution, exhibiting a 
delay-induced reduction of the spontaneous spiking activity at characteristic frequencies. Moreover, 
a specific frequency-locking mechanism is detected for the mean interspike interval. 

PACS numbers: 05.40.Ca, 87.18.Tt, 87.19.1m, 87.16.A- 
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I. INTRODUCTION 

Time-delayed feedback is a common mechanism rele- 
vant in many biological system including excitable gene 
regulatory circuits [l[ and human balance [2j . It has been 
investigated not only in biological systems, but also in 
a wider area: This includes but is not limited to semi- 
conductor superlattices [H, chemical oscillators (CO ox- 
idation on platinum) [1] and photosensitive Oregonator 
models 5]. These types of systems have been formu- 
lated mathematically in terms of recurrent models @, 0] ■ 
Recurrent connections introduce delayed feedback loops 
which may dramatically change the dynamic behavior of 
the system. It is known that delayed feedback presents 
an efficient method to control chaos or turbulence via 
stabilizing the unstable periodic orbits (UPOs) embed- 
ded in the chaotic attractor [§[ . Feedback can be used to 
stabilize periodic orbits or to control the coherence 
resonance [!, [H, [T3| ; the latter has also been studied ex- 
perimentally for bistable system with delayed feedback 

Over the past decades, neurobiologists found out that 
axons propagate the neuron's electrical signal to other 
neurons, and may sometimes feedback to the same neu- 
ron's dendrites [12h15{ . These auto-synapses which es- 
tablish a time-delayed feedback mechanism on a cellular 
level, are called autapses and were described by Van der 
Loos and Glaser in 1972 [16]. Since the discovery of such 
autapse in pyramidal cells in the cerebral neocortex, they 
were found in about 80% of all analyzed neurons includ- 
ing neurons of the human brain |17| . However, what 
functional significance they offer in the neural systems 
is still not fully understood, as the autapses exhibit a 
broad range of delay time scales: these range from a few 
milliseconds to tenths of milliseconds Q- 

It was established that delayed feedback can induce 
bursting [l(| ■ Its role for coding and processing of infor- 
mation in the brain has been evidenced experimentally 



for a variety of neural systems. However, the mechanisms 
leading to bursting are still open to debate; in particular 
the presence of noise is expected to play a key role. 

Within this work we aim to contribute to this objec- 
tive by considering the influence of intrinsic noise. It has 
been demonstrated, that noise leads to nontrivial effects 
in neuronal dynamics [l8[. Som e typic al examples are 
stochastic resonance phenomena [19M22j j . coherence res- 
onance [23 26] and dynamical synchronization phenom- 
ena [27l - l30j ] . In our situation there occurs an intrinsic 
source of noise which is due to the stochastic gating of 
ion channels; i.e. so-called channel noise. The latter is 
inherently coupled to the properties of the axonal cell 
membrane and, a priori, cannot be neglected |3l| . Inter- 
estingly, this intrinsic noise effects the neuronal signaling 
at different levels: It was shown that it may control the 
occurrence of spontaneous action potentials, can improve 
the output of signal quality [2l[ and may also account for 
the reliability of propagation [321 ] , to name but a few. 

The manuscript is organized as follows: in Sec. [H] we 
present the biophysical model: It is given in terms of 
a stochastic Hodgkin-Huxley model containing a time- 
delayed feedback current. In Sec. lIIII we discuss the repet- 
itive firing that occurs for the deterministic dynamics 
(i.e. in absence of noise) of this Hodgkin-Huxley model 
with delayed feedback. The stochastic dynamics are ad- 
dressed in Sect. IIV1 Finally, we present our conclusions 
with Sec. El 



II. MODEL SET UP 

In 1952, Hodgkin and Huxley proposed an archetyp- 
ical model for cell excitability and signal transmission 
along the axon of a neuronal cell [33|. They postulated 
a set of continuous parallel pathways for the passage of 
the ionic and capacitive currents to represent the elec- 
trical properties of the neuronal cell membrane. The 
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FIG. 1: Sketch of an autapse: Sketch of a neuronal cell ex- 
hibiting a self delayed feedback mechanism. 



Hodgkin-Huxley model is widely regarded as a milestone 
achievement in biophysics and especially in electrophys- 
iology. Its applicability was extended to more complex 
systems beyond the originally studied excitability of the 
squid giant axon. Experimental evidence of the chan- 
nel noise has resulted in stochastic generalizations of the 
Hodgkin-Huxley model by considering the stochastic dy- 
namics of the ion channel gating [2~3| . These generaliza- 
tions allow for microscopic modeling of the occurrence 
of spontaneous spiking and the phenomenon of biological 
Stochastic Resonance [20|. In this context, it was shown 
recently, that channel noise improves the reliability for 
signal transmission along neuronal axons (32j. 

We mathematically model the delayed feedback in- 
duced by an autapse (sketched in Fig. Q]) by a feedback 
mechanism within a Hodgkin-Huxley model of the type 
put forward by Pyragas 35] . By an additional pathway in 
the dynamics of the membrane potential we consider the 
autapse leading to a delayed stimulus I T (t) . Particularly, 
we consider an excitable dynamics for the membrane po- 
tential V(t), reading: 



c av(t)_ = n 4 (f)5 ™* {v{t) _ Vk) _ gL {v{t) Vl) 



dt 



- m 3 (t)h(t)g^ (V(t) - K Na ) + J r (t) 
wherein the autaptic delayed stimulus 

I r (t) =e[V(t-r)-V(t)} , 



(1) 



(2) 



is proportional to the difference of the membrane poten- 
tial at time t and that at an earlier time t — r. Here, 
e corresponds to the coupling strength of the autapse 
mechanism and the finite delay time r refers to the spe- 
cific time delay caused by the autapse which is due to 
a finite signal propagation speed. The delay time r is 
representing the elapsed time associated with the axonal 
propagation prior to the signal recurring onto the neu- 
ron; V(t — t) is the membrane potential at the earlier 
time t-r. 

The autaptic delayed stimulus in Eq.[5]results in an ex- 
citatory coupling mechanism in which spiking of a cell at 



an earlier time t— t favors the initiation of a spiking event 
of the same cell at time t. We note, that the ansatz for the 
delayed self-stimulus corresponds to electrotonic interac- 
tion, i.e. we consider an idealized situation wherein the 
autaptic delayed stimulus is proportional to the differ- 
ence of the pre- and post-synaptic membrane potentials. 
Moreover, as autapses are typically formed by chemical 
synapses, our modeling simulates the complex biophys- 
ical temporal evolution of the synaptic conductance by 
invoking a constant coupling strength e. 

In Eq. |T]), C denotes the membrane capacitance, V(t) 
the time-dependent membrane potential, Vl the leak- 
age potential, Vk and V^a the reversal potentials for 
the potassium and sodium currents. The leakage con- 
ductance is given by gK, the potassium and sodium 
maximum conductances read g£? ax and g™*, respec- 
tively. The functions m(t), n(t) and h(t) in Eq. ([1]) 
denote the so-called stochastic gating variables (see be- 
low), describing the mean fraction of open gates of the 
sodium and potassium channels at time t. The typi- 
cal values of parameters in our numerical simulation are 
those of the original Hodgkin-Huxley model [HI; i.e., 
(7=1 ^F/cm 2 , g™ ax = 120mS/cm 2 , g™ x = 36mS/cm 2 , 
g h = 0.3mS/cm 2 , V K = -77mV, V Na = 50mV and 
V L = 54.4mV. 

The stochastic dynamics of the gating variables m(t), 
n(t) and h(t) depend on the voltage-dependent opening 
and closing rates oti{V) and Pi{V) (i = m,h,n), which 
read 1331] : 



ot m {V) = 

Pm(V) = 

a h (V) = 



0.1(V + 40) 



a n {V) 



l-exp{-(F + 40)/10} ' 
4 exp{- (F + 65)/18} , 
0.07 exp{-(y + 65)/20} , 
1 

l + exp{-(T/ + 35)/10} ' 

0.01(7 + 55) 
l-cxp{-(^ + 55)/10}' 
0.125 exp{-(^ + 65)/80} 



(3a) 

(3b) 
(3c) 

(3d) 

(3e) 
(3f) 



The resulting stochastic dynamics is then modeled by 
a Fokkcr-Plank-type dynamics for the individual gates. 
Specifically, this dynamics emerges within a large sys- 
tem size-expansion of the underlying Markovian master- 
equation dynamics for the dynamics of the number of 
open gates [HI, HH . This so obtained Langevin-dynamics 
is then interpreted in the Ito sense, reading explicitly: 



di 
dt 



Oi{V)(l -i)-0i(V)i + Si(t), 



(4) 



with i — m,h,n. Here, £j(t) (i — m,h,n) denote in- 
dependent Gaussian white noise sources with vanishing 
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FIG. 2: The phase diagram e— r for repetitive firing within the 
Hodgkin Huxley model containing the delay coupling strength 
e is depicted vs. the delay time r. The two black solid lines 
indicate the boundaries of critical coupling strength e versus 
the delay time r at which the behavior changes from the rest 
state to a repetitive firing. The gray areas give the parameter 
regime for which repetitive firing is observed. 



mean and auto-correlation functions: 

(U(t)UV)) = (1 ~ m) " m+W/3m ^ - t% (5a) 

mw)) = {1 ~ h) " h+hph s(t-n ( 5b) 

(€n(*)€n(*0> = (l ~ n ^ + n/ V - f )■ (SC) 

Note, that the noise strengths are determined by 
the number of potassium and sodium channels Nk and 
A^Na- Throughout this work we assume constant ion 
channel densities; i.e. following the original Hodgkin- 
Huxley model we use 18 for potassium and 60 for the 
sodium channels per /im 2 [33J. An integration step 
At = 0.001ms was used in the simulations and for the 
generation of the Gaussian distributed random numbers, 
the Box-Muller algorithm [33] was used. 



III. DETERMINISTIC DYNAMICS 



delay time parameter r, there emerges now a stable, os- 
cillatory solution, presenting repetitive (tonic) spiking. 
To investigate this repetitive spiking behavior in more 
detail, we systematically varied the two parameters (e 
and r) of the delay coupling scheme. In doing so, we de- 
termined the first occurrence of repetitive firing after a 
spike was created. The resulting phase diagram with the 
positive and negative critical coupling strength e+ (r) and 
£~(t), respectively, is depicted in Fig. [2j We point out 
the resulting asymmetry between positive and negative 
coupling e. 

For sub-threshold coupling strength, i.e. e~(r) < e < 
e+(r) (in the white regions of Fig. [5]), excitations die 
out and the system relaxes to the rest state. For supra- 
threshold coupling, i.e. e > e+(r) or e < e~(r) a repeti- 
tive firing is observed. If the delay time becomes shorter 
than the refractory time interval, being ca. 12 ms, which 
is the time the system needs for spiking and returning 
to the rest state, the modulus of the critical coupling 
strength of the bifurcation increases drastically. This is 
due to the fact, that in the so-called undershoot phase 
(i.e. the part of the course of the action potential when 
the membrane potential is smaller than the membrane 
potential at rest) the neuron is rather insensitive to any 
stimuli. Hence, a stronger coupling strength is needed 
in order to excite the system from the undershoot phase 
and to obtain repetitive firing. In case of larger delay 
the value of the critical coupling saturates. For longer 
delay times more than one spike may fit into time in- 
terval given by r. Consequently, doublets, triplets and 
multiplets may appear. 

Below, we shall present examples where we have chosen 
r = 35 ms and used a positive coupling strength. The 
critical value for repetitive firing corresponds to e+ = 
0.059 mS/cm 2 . For larger coupling e > e+(r) an action 
potential induces repetitive firing. In section HVl we refer 
to this regime as supra-threshold delay coupling regime. 
In case of lower coupling strength; i.e. e < e+(r) the state 
of a repetitive firing is excitable and noise is needed to 
support the repetition of an action potential. In the next 
section HVl this regime will be addressed as sub-threshold 
delay coupling regime. Remarkably, for negative coupling 
e < one finds an overall similar similar behavior that 
differs, however, quantitatively (not shown). 



In the absence of both, the delayed feedback, i.e. e = 0, 
and zero noise strength, with the latter being formally 
achieved by N^ a — ► oo and Nk — > oo, the original 
Hodgkin-Huxley dynamics are recovered. These exhibit a 
single stable fixed point, which is the rest state with the 
rest potential Vr es t ~ — 65mV. Any temporary distur- 
bance, e.g. caused by an applied external current stimu- 
lus, decays and the system relaxes back to the rest state. 
Note that upon a proper choice of the initial condition, 
a single action potential may occur, before the system 
relaxes to the rest state (38j . 

In presence of the autaptic delay coupling, the fixed 
point solution remains stable, but upon increasing the 



IV. STOCHASTIC DYNAMICS OF THE 
HODGKIN-HUXLEY MODEL WITH DELAYED 
FEEDBACK 

In Fig. [3] the simulated membrane potential V(t) is 
shown for an exemplarily chosen noise level and for two 
cases, namely in panel (a) without delay coupling (e = 0) 
and in panel (b) with a finite delay coupling (e > 0). In 
absence of delay coupling the occurrence of spontaneous, 
i.e. noise induced actions potentials occur irregular while 
a characteristic bursting pattern can be detected in pres- 
ence of finite delay-coupling e. 
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FIG. 3: Simulated spike-trains: (a) Membrane potential V(t) 
from the stochastic Hodgkin-Huxley model in absence of feed- 
back, i.e. with e = versus time t; (b) V(t) versus time t in 
presence of finite feedback of strength e = 0.07mS/cm 2 ex- 
hibiting repetitive spiking for a chosen delay time r = 35 ms. 
The number of potassium ion channels is set to Nk = 300 
and those of the sodium ion channels to A^Na = 1000. The 
spontaneous, i.e. noise induced, action potentials exhibit a 
bursting-like behavior in case of the delay coupling (b). In 
(a) we indicated with the double-arrow line a typical occur- 
ring noise-induced interspike interval TV 



In order to explain the spiking dynamics quantita- 
tively, we determined the time intervals between two 
spike events Tj, see in Fig. EJa), as obtained from sim- 
ulations of the stochastic Hodgkin-Huxley model with 
delayed feedback coupling, cf. Eq. (|TJ). Note that in 
order to detect the occurrences of action potentials from 
the simulated membrane potential dynamics V(t), we did 
define a specific threshold barrier for detection. In par- 
ticular, whenever the membrane potential V(t) exceeds 
the value of mV, the occurrence of an action potential 
is assigned. In fact, the so determined spike occurrences 
depend only weakly on the actual choice of the detection 
barrier (2lj . 

These observed stochastically emerging interspike in- 
tervals Ti the render the (normalized) interspike interval 



histograms (ISIH) and the mean interspike interval (T) 

1 N 

CO = ^£r i5 (6) 
1=1 

where TV denotes the number of spikes obtained in the 
individual numerical simulation. The bin width for the 
histograms has been set at 0.2 ms. Both these quantities 
form the basis of our analysis of spike trains. 

We begin this study with the dynamics of zero delay 
coupling for the neuron, i.e. with e = 0. The distribution 
of the interspike intervals at a fixed channel noise inten- 
sity is depicted in Fig. @] For large time intervals, the 
distribution decays exponentially. Small time intervals 
are suppressed because of the neuron's refractory state. 
For e = the ISIH exhibits one broad maximum located 
around the internal time scale TTnt • This time scale sys- 
tematically diminishes with increasing noise level [IH l2l| . 
The minimal time scale for undistorted spikes is around 
the refractory time of ~ 12ms. 

In presence of noise and finite feedback coupling, the 
neuron still can return to its fix-point with nonvanishing 
probability. This is why the peak at the intrinsic time 
scale and a broad exponential decay is still detected in 
the ISIH in Fig. 0J However, the ISIHs become multi- 
modal in case of finite, nonvanishing feedback. Several 
maxima occur at larger time intervals and, in between, 
a number of deep dips are observed. Such forms of the 
distributions of the ISIH are indicators of an enhanced 
coherence [TH, [H| . Whereas in the sub-threshold regime 
(dashed line) the additional extrema merely show up, 
they become increasingly dominant in the case of supra- 
threshold coupling (note the logarithmic scale of the or- 
dinate in Fig. HJ. 

A. Sub-threshold delay coupling: e < ej (r) 

In the case of sub-threshold coupling the repetitive os- 
cillatory spiking is noise supported, as a non- vanishing 
delayed stimulus I T (t) shifts the stable fix-point towards 
the threshold for excitation. Consequently, the gener- 
ation of spikes becomes more likely. However, the dy- 
namics is still excitable; i.e. a threshold still exists and 
noise is needed to induce its excitation. Consequently, a 
sub-threshold delayed stimulus favors some characteristic 
interspike interval time T T - As far as the system remains 
excitable, an activation time T ac t after the delayed stim- 
ulus did set in is necessary in order to create the next 
noise-induced spike. Therefore, the time T T is not equal 
to the delay time r, but equals the sum of the delay time 
and an activation time 7Ict> being around 2 ms in this 
considered case. 

In addition, there occurs a suppression of certain time 
scales when the delay coupling is applied to the neuronal 
dynamics, cf. Fig. For positive e the system favors 
frequencies associated to the time T T and suppresses fre- 
quencies slightly larger than those. The feedback hinders 
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FIG. 4: (Color online) Interspike interval histograms (ISIH) 
versus the interspike interval T of the stochastic Hodgkin- 
Huxley model with a delay coupling mechanism at different 
coupling strengths e, as indicated in the figure in units of 
mS/cm 2 . The bin width for T has been chosen throughout 
this work at 0.2 ms. The chosen number of potassium ion 
channels is Nk = 150, while the number of sodium ion chan- 
nels is set at A^a = 500. The delay time r is set at 35 ms. 



FIG. 6: (Color online) Interspike interval histogram (ISIH) for 
different values of the supra-threshold coupling strength e > 
et(r): The distribution of the interspike intervals is depicted 
for fixed delay time r = 35 ms, the numbers of ion channels 
are Nk = 150 and A^Na = 500. The values of the coupling 
strengths indicated in the plot are in units of mS/cm 2 . 
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FIG. 5: (Color online) Interspike interval histogram (ISIH) vs. 
interspike duration T for different noise levels corresponding 
to different number of sodium and potassium channels as in- 
dicated. The chosen delay time is : r = 35 ms; the coupling 
strength is set at e = 0.03mS/cm 2 < ej(r) and Nn& = -f Wk. 



the generation of spikes by noise if the feedback current 
I T {t) approaches values around the polarized refractory 
period. 

In the sub-threshold regime the intrinsic time scale 7I n t 
and correspondingly the height and the position of the 
associated peak distinctly depends on the noise level. In 
Fig. [5] we depict results from simulations for different 
noise levels, i.e. different numbers of embedded ion chan- 
nels. With increasing noise level the generation of spon- 
taneous spikes by noise gains influence substantially. As 
a result, the intrinsic time scale becomes changed and 
the larger time intervals reflect the broad, exponentially 
decaying ISIH. The former sharply peaked maximum at 
T T switches towards a noise supported broader peak at 



smaller times, and the maximum shifts with growing 
noise towards the minimal time for a interspike inter- 
val located at the refractory time. Interspike intervals 
around this minimal time are promoted mainly by the 
presence of noise. They do not necessarily relate to the 
delay mechanism. In contrast, the time scale T T induced 
by the delayed feedback is only slightly affected by the 
noise level. This is a somewhat striking feature recall- 
ing the fact that in the case of an excitable dynamics 
these spikes in fact require the presence of noise to be- 
come activated. The robustness of the positions of the 
time scales as well as their separation from the noisy 
time scale in terms of the resulting bimodal shape of the 
ISIH thus indicates a noisy synchronization phenomenon 

[IS 13, 53- 



B. Supra-threshold delay coupling: e > e+(r) 

We next consider the case with coupling strengths be- 
yond the critical strength, i.e., e > e+(r). Here, the 
spikes repeat deterministically. However, noise-induced 
skipping of spikes [39j, |40j leads to a transition from the 
oscillatory state with repetitive firing to the excitable 
state. In turn, the back-transitions from the excitable to 
the oscillatory state are caused by noise-induced sponta- 
neous spikes, cf. Fig. [3J The spikes repeated determinis- 
tically lead to sharp peaks in the ISIH, while the noise- 
induced back-transitions from excitable to the oscillatory 
state result in a broad distribution with exponentially 
decaying tail. In Fig. [SJ the resulting ISIH is depicted 
for various supra-threshold coupling strengths. Upon the 
chosen parameter values for the noise level, i.e. the num- 
ber of sodium and potassium ion channels, and the cou- 
pling parameters e and r, the ISIH therefore exhibits 
sharp peaks and more or less pronounced broad back- 
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FIG. 7: (Color online) Average interspike interval (T) as 
function of the delay time r: The mean interspike inter- 
val (T), cf. Eq. ([6]), is plotted versus the delay time r for 
the following parameters: the used ion channel numbers are 
Nk = 150, JVisra = 500 and the coupling strengths are set 
at e = 0.1mS/cm 2 (black dashed line), e = 0.2mS/cm 2 (red 
dotted line), and e = 0.4mS/cm 2 (blue solid line). 



ground: For stronger coupling or weaker noise the broad 
background diminishes. Contrarily, the peak height at T T 
shows a strong dependence on the noise level, indicating 
the competitive interplay between channel noise and the 
delayed feedback mechanism: With increasing noise level 
the height of the peak at T T decreases, while the intrinsic 
time scale acquires increasing influence (depicted for the 
sub-threshold coupling regime in Fig. [5]). Note that we 
have selected the delay time (r = 35 ms) so that the sep- 
aration between the broad background which is peaked 
at the intrinsic time scale 7Tnt and the delayed induced 
sharp peaks at T T and multiples thereof is educible and 
visible with the ISIHs, cf. Fig. [5] and Fig.[BJ 

Strikingly, with increasing coupling strength the mul- 
timodal structure collapses to a unimodal one, cf Fig. [5] 
The most probable interspike interval now centers at one 
value. Due to this supra-threshold driving, each sponta- 
neous spiking event is repeated periodically and even the 
noise dominated scale now collapses towards a narrow 
peak. Note also that the distribution gains in sharpness 
as the coupling strength increases, cf. Fig. [6] 

Let us focus on the dependence of the interspike inter- 
vals on the delay time r. For this purpose we consider 
the mean interspike interval, cf. Eq. ((6]). This quantity 
inherits the characteristic dependence on r since the ISIH 
shape is unimodal. Larger deviations from the mean in- 
terspike interval are found during transitions between dif- 
ferent synchronized states where the histogram becomes 
multimodal. 

In Fig. [7j the mean interspike interval (T) is depicted 
for a fixed channel noise strength. If the delay coupling 
mechanism does not dominate the spiking, i.e. for small 
supra-threshold coupling strengths, the mean interspike 
interval exhibits a smooth dependence on the delay time. 
However for larger supra-threshold driving strengths, for 




FIG. 8: The inverse of the number of spikes that fit into the 
delay-time interval k is plotted as function of the delay time 
r for TVk = 150, -/Vk = 500 and a delay coupling strength of 
e = 0.4mS/cm 2 . 



which the ISIH exhibits a unimodal structure consisting 
of a sharp peak, the mean interspike time varies with the 
delay time r in an almost piecewise linear fashion, dis- 
playing sharp, triangle-like textures , cf. Fig. [7J At these 
sharp peak locations, the number of spikes that fit ac- 
cording to the intrinsic time into a full time length given 
by the delay time, just increases by unity with increasing 
delay time r. The mean interspike interval is henceforth 
proportional to the ratio of the delay time and the num- 
ber of spikes n fitting into this very delay-time interval; 
yielding 



(T) oc t / n = 



(7) 



In Fig. [5] the behavior of k vs. the delay time r is de- 
picted. At multiples of the noise-dependent intrinsic time 
scale, such characteristic steps do occur indeed. 

We find that stronger noise intensity mimics a decrease 
of the coupling strength. Nearby those typical steps noise 
is able to induce newly generated spikes at the end of the 
intervals reducing thereby the ratio k. Alternatively, the 
role of noise decreases the number of spikes towards the 
left boundary of a new synchronized region. Therefore 
the corresponding results at increased noise strength at- 
tain a form that is similar to the case of lower delay 
coupling strength. 

A comparable scenario has been reported in other sys- 
tems such as for a noisy Van der Pol dynamics, see in 
Refs. [42T - |44T | near the Hopf bifurcation, for semiconduc- 
tor superlattices Q and stochastic excitable dynamics 
@. We find that noise can induce oscillations with a 
well-defined time scale, and the simultaneous application 
of a delay of duration r is able to stabilize this orbit. 



C. Phase oscillator modeling 

We have seen that the neuron behaves almost like an 
oscillatory system whenever the delay coupling rules the 



7 



dynamics. This brings to mind a description in terms of 
a phase oscillator model which we develop next. Let us 
assign an intrinsic frequency by oj — 27r/7i n t- As pointed 
out above, the autaptic feedback leads to a stabilization 
of the internal rhythmicity as determined by the model 
parameters. This in turn results in synchronized pat- 
terns, cf. Fig. m 

Therefore we investigate the question whether the ob- 
served synchronization in this studied neuronal model 
with a utap se can be characterized by a Kuramoto-type 
model 45] of a phase oscillator including delayed feed- 
back? Put differently, upon neglecting the details of the 
shape of the spike trains the essential information can be 
reduced to the phase of the oscillatory neuron. 

Let us next compare our findings with the results ob- 
tained from a pure phase dynamics. Accordingly, the 
proposed dynamics for the underlying phase is modeled 



4> = oj - A sin [4>(t) - <j){t - t)] 



(8) 



where oj denotes the angular intrinsic frequency of the 
oscillator and A is the delay-coupling strength. Searching 
for solutions with fixed frequency O, we make the Ansatz: 

(t>(t)=flt. (9) 

The resulting, self-consistent equation reads [46T - I48I ] : 

n = oj - A sin [£It] , (10) 

which was solved numerically. Comparing this phase os- 
cillator dynamics with our oscillatory stochastic neuron 
dynamics exhibiting an autapse, we make the following 
parameter identification, i.e., Q — 2tt/(T) and obtain 
from an expansion in small delay time r the relation 
A = l/7Int- In Fig.[5]we present the comparison between 
the Hodgin-Huxley modeling and the the phase oscilla- 
tor modeling: Even though the quantitative behavior is 
different, a good qualitative agreement is detected. The 
positions of the spikes where the autapse model jumps be- 
tween bursts of different number of multiplets are nicely 
reproduced by this simplified phase oscillator modeling. 
A further improvement of the modeling would require a 
more optimal determination of the r-dependent parame- 
ters of this Kuramoto model with feedback. 



V. CONCLUSION 

With this work we have investigated the effect of de- 
layed feedback on the dynamics of a stochastic Hodgkin- 
Huxley neuron. Using the original Hodgkin-Huxley 
parameters, we simulated numerically the stochastic 
Hodgkin-Huxley model taking into account the effects 
of intrinsic channel noise. A Pyragas-like delayed feed- 
back mechanism is employed to model the autapse phe- 
nomena, in which a neuron's dendrite back-couples to 
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FIG. 9: (Color online) Comparison of the stochastic, neuronal 
dynamics Hodgkin-Huxley model including feedback, Eqs. fT] 
OH [4] [5]) with the dynamics of the feedback assisted Kuramoto- 
type phase oscillator model, Eq. (|SJ. The mean interspike 
interval (red dashed lines) is the result of the numerical inte- 
gration of the stochastic Hodgkin-Huxley model with delayed 
feedback; there is quantitative good agreement between pe- 
riod of oscillations as obtained from a phase oscillator dynam- 
ics of the Kuramoto-type delayed feedback model (solid black 
line). Regions with multiple solutions in the phase model 
correspond to discontinuous jumps between different types of 
multiplets in the Hodgkin-Huxley model. The parameters for 
the neuronal modeling are chosen as Nk = 150, JVnb — 500 
and e = 0.4mS/cm 2 ; those of the phase oscillator are given 
by oscillator frequency ui — 2-7r/15.5 ms and feedback strength 
A = 1/15.5 ms. 



itself. The two basic parameters used in our study are 
the strength of the autapse-coupling e and the time-delay 
r resulting from the finite length of the self-connecting 
dendrite. 

We have found that the neuron dynamics exhibits in- 
triguing new times scale that stem from the autaptic con- 
nection. Since the rest state of the neuron is always sta- 
ble, noise or an initial spike is necessary to create activity, 
i.e. the spiking dynamics. For a small number of Na and 
K channels, the noise becomes sizable and the excitory 
dynamics remains practically unaffected by the delay. 
In contrast, smaller noise levels and stronger coupling 
strengths induce a different synchronization phenomena 
between the delay time and the intrinsic (also noise de- 
pendent) time scales. The delay time and these intrinsic 
time scales determine how many spikes will be created 
and become subsequently locked during one delay epoch. 
We further underline that our exemplary study can be 
mimicked qualitatively in terms of a reduced description 
given by a Kuramoto phase dynamics with built-in feed- 
back. 

We have shown, that the delayed feedback mechanism 
serves as a control option for adjusting the peaked dis- 
tribution of interspike intervals, being of importance for 
memory storage [49| and stimulus-locked short-term dy- 
namics in neuronal systems [50| . One may therefore spec- 
ulate whether nature adopted the autapse phenomena for 
frequencies-filtering in presence of unavoidable intrinsic 
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channel noise. 
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